An extended Hubbard model with ring exchange: a route to a non-Abelian topological phase 
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We propose an extended Hubbard model on a 2D Kagome lattice with an additional ring exchange term. The 
particles can be either bosons or spinless fermions. At a special filling fraction of 1 /6 the model is analyzed 
in the lowest non-vanishing order of perturbation theory. Such an "undoped" model is closely related to the 
Quantum Dimer Model. We show how to arrive at an exactly soluble point whose ground state is the "d- 
isotopy" transition point into a stable phase with a certain type of non-Abelian topological order. Near the 
"special" values, d = 2 cos 7r/(fc + 2), this topological phase has anyonic excitations closely related to SU(2) 
Chern-Simons theory at level k. 



Since the discovery of the fractional quantum Hall effect 
in 1982 d, topological phases of electrons have been a 
subject of great interest. Many Abelian topological phases 
have been discovered in the context of the quantum Hall 
regime I2I1. More recently, high-temperature superconduc- 
tivity fsv?, 9, ^, "t", "s", ^ and other complex materials have 
provided the impetus for further theoretical studies of and ex- 
perimental searches for Abelian topological phases. The types 
of microscopic models admitting such phases are now better 
understood |10, ll^Jl- 

Much less is known about non-abelian topological phases, 
apart from some tantalizing hints that the quantum Hall 
plateau observed at v = 5/2 might correspond to such a non- 
abelian phase II5L II6II . However, non-abelian topo- 
logical states, if created and controlled, would open the door 
to scalable quantum computation |17, 18]. As a first step, 
the study of a class of topological field theories has been re- 
duced to combinatorial manipulations of loops on a surface 
1 19, 20]. A virtue of this formulation is that it exposes a strat- 
egy for constructing microscopic physical models which ad- 
mit the corresponding phases; since Hilbert space is reduced 
to a set of pictorial rules, the models should impose these rules 
as energetically favorable conditions satisfied by the ground 
state. In this paper, we show how this approach can be imple- 
mented. 

We propose a microscopic model which has the following 
properties: (a) it is an extension of the Hubbard model and, 
therefore, is quasi-realistic, (b) it is soluble, (c) for certain 
model parameters, it is perched at a transition point f2lll into 
a non-Abelian topological phase relevant to quantum com- 
putation. By quasi-realistic, we mean that the model has 
short-ranged interactions and hopping, so it is possible that 
the Hamiltonian of a real material could be viewed as a small 
perturbation of the Hamiltonian of this paper. Optical lattices 
12211 . quantum dot or Josephson junction arrays 12311 might be 
designed with Hamiltonians in this general class, and these 
may also be promising avenues for realizing our model. 

The non-Abelian topological phases refeiTed to in the above 
paragraph are related to the doubled SU(2)fe Chern-Simons 
theories described in ]20, 24]. These phases are characterized 
by (fc + l)^-fold ground state degeneracy on the torus and 
should be viewed as a natural family containing the topolog- 



ical (deconfined) phase of Z2 gauge theory as its initial ele- 
ment, k = 1. For k > 2 the excitations are non-Abelian. For 
k = 3 and fc > 5 the excitations are computationally universal 
Here, we describe the conditions which a microscopic 
model should satisfy to be in such a topological phase. It is 
useful to think of such a microscopic model as a lattice reg- 
ularization of a continuum model whose low energy Hilbert 
space may be described as a quantum loop gas. More pre- 
cisely, a state is defined as a collection of non-intersecting 
loops, as discussed in ]21, 24, 26]. A Hamiltonian acting on 
such state can do the following: (i) the loops can be continu- 
ously deformed - we will call this "move" an isotopy move; 
(ii) a small loop can be created or annihilated - the combined 
effect of this move and the isotopy move has been dubbed 'd- 
isotopy 'f20l|2ll|2l; (iii) finally, when exactly fc + 1 strands 
come together in some local neighborhood, the Hamiltonian 
can cut them and reconnect the resulting "loose ends" pair- 
wise so that the newly-formed loops are still non-intersecting. 
More specifically, in order for this model to be in a topological 
phase, the ground state of this Hamiltonian should be a super- 
position of all such pictures with the additional requirements 
that (i) if two pictures can be continuously deformed into each 
other, they enter the ground state superposition with the same 
weight; (ii) the amplitude of a picture with an additional loop 
is d times that of a picture without such loop; (iii) this super- 
position is annihilated by the application of the Jones- Wenzl 
(JW) projector that acts locally by reconnecting fc + 1 strands. 
Readers interested in the details are referred to ]20, 24, 2^ 
and references therein. In this paper, we focus on the first two 
conditions, which place the system at a transition point into 
the desired phase(s) ] 2 1 ] ; our purpose is to construct a Hamil- 
tonian which enforces d-isotopy for its ground state(s). 



Our proposed model is defined on the Kagome lattice 
shown in Fig. ^ The sites of the lattice are not completely 
equivalent, in particular we choose two special sublattices - TZ 
(red) and Q (green) whose significance will be discussed later. 
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FIG. 1 : Solid dots and dashed lines represent sites and bonds of the 
Kagome lattice K, with the special sublattices IZ (red) and Q (green). 
Solid lines define the surrounding triangular lattice. 

The Hamiltonian is given by: 

H = ^^^jLini+Uo'^n1+U ninj+ ^ VijUiUj 

i i {iJ)eO (i,i)eix],^0 

-E*^j(c1cj + 4c,) + tr. (1) 

Here rii = cjc^ is the occupation number on site i, /ij is the 
coiTesponding chemical potential. Uq is the usual onsite Hub- 
bard energy Uq (clearly superfluous for spinless fermions). U 
is a (positive) Coulomb penalty for having two particles on the 
same hexagon while Vij represent a penalty for two particles 
occupying the opposite corners of "bow ties" (in other words, 
being next-nearest neighbors on one of the straight lines). We 
allow for the possibility of inhomogeneity so not all Vij are 
assumed equal. Specifically, define wjjj, = Vij where a is the 
color of site (i), b is the color of (j), and c is the color of the 
site between them. In the lattice in Fig. ^ we have, possibly 
distinct, v^^, v^f^, v^^, v'f.^^, and v'^.g, where r € TZ, g € G, 
and b E B = /C\(7?. U G)- Uj is the usual nearest-neighbor 
tunnelling amplitude which is also assumed to depend only on 
the color of the environment: tij = where c now refers to 
the color of the third site in a triangle. Finally, we include Tr 
- a four-particle ring exchange term whose exact form will be 
specified later. Tr is added to the Hamiltonian on an ad hoc 
basis to allow correlated multi-particle hops. Ring exchange 
terms can be justified semiclassically iEtIi and they do indeed 
appear in in such physical systems as spin systems |28, 29], 
solid '^He llsoll and Wigner crystals |31]. Of course, small ring 
terms can arise perturbatively along the lines of |29], e.g. a 
four-particle move occurs at order 4. 

The onsite Hubbard energy Uq is considered to be the 
biggest energy in the problem, and we shall set it to infinity, 
thereby restricting our attention to the low-energy manifold 
with sites either unoccupied or singly-occupied. The rest of 
the energies satisfy the following relations: U 2> Uj, Vij, fii; 
we shall be more specific about relations between various 
tij's, Vij's and fii's later. 



The "undoped" system corresponds to the filling fraction 
1/6 (i.e. Np = J2i "-i = ^/ 6, where N is the number of sites 
in the lattice). The lowest-energy band then consists of con- 
figurations in which there is exactly one particle per hexagon, 
hence all ?7-terms are set to zero. These states are easier to 
visualize if we consider a triangular lattice T whose sites co- 
incide with the centers of hexagons of /C. (K, is a surrounding 
lattice for T.) Then a particle on K, is represented by a dimer 
on T connecting the centers of two adjacent hexagons of /C. 
The condition of one particle per hexagon translates into the 
requirement that no dimers share a site. In the 1/6-filled case 
this low-energy manifold coincides with the set of all dimer 
coverings (perfect matchings) of T. The "red" bonds of T 
(the ones corresponding to the sites of sublattice TZ) them- 
selves form one such dimer covering, a so-called "staggered 
configuration". This particular covering is special: it contains 
no "flippable plaquettes", or rhombi with two opposing sides 
occupied by dimers (see Fig.Q. 

So henceforth particles live on bonds of the triangular lat- 
tice (Fig.Q and are represented as dimers L341 . In particular, a 
particle hop corresponds to a dimer "pivoting" by 60° around 
one of its endpoints, Vij = is now a potential energy of 
two parallel dimers on two opposite sides of a rhombus (with c 
being the color of its short diagonal). It is clear that our model 
is in the same family as the quantum dimer model fT], which 
has recently been shown to have an Abelian topological phase 
on the triangular lattice lHoll which, corresponds to fc = 1, or 
d = 1. Here, we show how other values of k can be obtained. 

The goal now is to derive the effective Hamiltonian acting 
on this low-energy manifold represented by all possible dimer 
coverings of T. Our analysis is perturbative in t/U =: e. The 
initial, unperturbed ground state manifold for Uq — oo, U 
large and positive, all tij , Vij — and all fii equal is spanned 
by the dimerizations V of the triangular lattice T. As we grad- 
ually turn on the Vs, u's, and Tr, we shall see what equations 
they should satisfy so that the effective Hamiltonian on T) has 
the desired d-isotopy space as its ground state(s). 

Since a single tunnelling event in V always leads to dimer 
"collisions" (two dimers sharing an endpoint) with energy 
penalty U, the lowest order at which the tunnelling processes 
contribute to the effective low-energy Hamiltonian is 2. At 
this order, the tunnelling term leads to two-dimer "plaquette 
flips" as well as renormalization of bare onsite potentials /i^'s 
due to dimers pivoting out of their positions and back. We 
always recompute bare potentials /i^'s to maintain equality up 
to errors 0{e^) among the renormalized Jli's. This freedom to 
engineer the chemical potential landscape to balance kinetic 
energy is essential to finding an exactly soluble point. 

Let us pause and discuss the connection between our quan- 
tum dimer model and a desked topological phase . It is an old 
idea (see e.g. js^l) to turn a dimerization (perfect matching) 
J' into a collection of loops by using a background dimeriza- 
tion TZ to form a 'transition graph' TZ^ J. It turns out that 
fixing TZ as in Fig.Q without small rhombi with two opposite 
sides red, as the preferred background dimerization we obtain 
the fewest equations in and also achieve ergodicity |33] under 
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a small set of moves. Unlike in the usual case, the background 
dimerization TZ is not merely a guide for the eyes, it is phys- 
ically distinguished: the chemical potentials and tunnelling 
amplitudes are different for bonds of different color 

Let us list here the elementary dimer moves that preserve 
the proper dimer covering condition: 

(i) Plaquette (rhombus) flip - this is a two-dimer move around 
a rhombus made of two lattice triangles. Depending on 
whether a "red" bond forms a side of such a rhombus, its di- 
agonal, or is not found there at all, the plaquettes are referred 
to, respectively, as type 1 (or 1'), 2, or 3 (see Fig.|2}. The dis- 
tinction between plaquettes of type 1 and 1 ' is purely direc- 
tional: diagonal bonds in plaquettes of type 1 are horizontal, 
for type 1' they are not. This distinction is necessary since 
our Hamiltonian breaks the rotational symmetry of a triangu- 
lar (or Kagome) lattice. 

(ii) Triangle move - this is a three-dimer move around a trian- 
gle made of four elementary triangles. One such "flippable" 
triangle is labelled 4 in Fig.|2l 

(iii) Bow tie move - this is a four-dimer move around a "bow 
tie" made of six elementary triangles. One such "flippable" 
bow tie is labelled 5 in Fig.|2] 




FIG. 2: Overlap of a dimer covering of T (shown in thick black) with 
the red covering corresponding to the special sublattice TZ. Shaded 
plaquettes correspond to various dimer moves described in the text. 
Green sublattice is not shown. 

To make each of the above moves possible, the actual 
dimers and unoccupied bonds should alternate around a cor- 
responding shape. Notice that for both triangle and bow tie 
moves we chose to depict the cases when the maximal possi- 
ble number of "red" bonds participate in their making (2 and 
4 respectively). Note that there are no alternating red/black 
rings of fewer than 8 lattice bonds (occupied by at most 4 
non-colliding dimers). Ring moves only occur when red and 
black dimers alternate; the triangle labelled 4 in Fig.|5]does 
not have a ring exchange term associated with it, but the bow 
tie labelled 5 does: 




Here is the correspondence between the previous smooth 
discussion and rhombus flips relating dimerizations J' of T. 



Our surface is now a planar domain with, possibly, periodic 
boundary conditions (a torus). A collection of loops is gener- 
ated by 7?. U (with the convention that the dimers of 7?. fl 
be consider as length 2 loops or bigons). What about isotopy? 
Move 2 certainly is an isotopy from TZ U to TZ U J^' but 
by itself, it does almost nothing. It is impossible to build up 
large moves from type 2 alone. So it is a peculiarity of the 
rhombus flips that we have no good analog of isotopy alone 
but instead go directly to d— isotopy. We should impose the 
following relations associated with moves of type 1 (1'): 

since we pass from zero to one loop in (|3}- Additionally, the 
ring exchange term (|3 annihilates the superposition of one 
and four loops; we therefore require that X = d^. 

Having stated our goal, we now derive the effective Hamil- 
tonian H : V ^ V on the span of dimerizations. The 
derivation is perturbative to the second order in e where e = 
tbb/U = t^b/U. Additionally, t%/U = cqe where Co IS a 
positive constant, while = o(e) and can be neglected in 
the second-order calculations. (In the absence of a magnetic 
field all fs can be made real and hence symmetric with re- 
spect to their lower indices. Also, we set U — 1 for no- 
tational convenience.) We account for all second-order pro- 
cesses, i.e. those processes that take us out of T> and then 
back to v. These amount to off-diagonal (hopping) processes 
- "plaquette" flips" or "rhombus moves" - as well as diago- 
nal ones (potential energy) in which a dimer pivots out and 
then back into its original position. The latter processes lead 
to renormalization of the bare onsite potentials Hi, which we 
have adjusted so that all renormalized potentials jli are equal 
up to corrections 0(e'^). The non-constant part of the effec- 
tive Hamiltonian comes from the former processes and can 

be written in the form: H = J2x,j Axj- where 

H-ij is a 2 X 2 matrix corresponding to a dimer move in the 
two-dimensional basis of dimer configurations connected by 
this move. Ajj = 1 if the dimerizations I, e I? are con- 
nected by an allowed move, Aij = otherwise. Therefore 
it suffices to specify these 2x2 matrices Hij for the off- 
diagonal processes. For moves of types ( 1 )-(3), they are given 
below: 
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We can now tune H to the "small loop" value d. We require 

^(1) = ^(1') ^ 1^ 

_i j as these moves change the num- 
ber of small loops by one (cf. Eq. (|3ji). Since a move of type 2 
is just an isotopy move, we require cx i) - Finally, 
H^^^ = provided fc > 1, since it represents a "surgery" on 
two strands not allowed for fc > 1. (For k = 1, on the other 
hand, H'^'^^ oc ( \ )■) level k ~ \ configurations which 
differ by such a surgery should have equal coefficients in any 
ground state vector ^I^ while at levels fc > 1 no such relation 
should be imposed. Thus, for fc > 1 the matrix relations ( l4al 
^} yield equations in the model parameters: 

Types (1)&(1') : v^, = vl, = 2dc„e' (5a) 
and i;^^ = t-^g = 2d-icoe2 (5b) 
Types (2)&(3) : = 2e^ and = (5c) 

We have already assumed that the Hamiltonian has a bare 
ring exchange term, Tr given by Eq. (|2j or, in matrix form, 

Tr = a ^ ^ where \ = according to the discussion 

after Eq. (|3}- Additionally, we would want the off-diagonal 
elements of Tr, — aA to be of order thus making sure that 
this ring exchange dominates all other ring exchanges that will 
appear in the higher orders of perturbation theory. Along with 
Eqs. Q, these conditions place our model at the soluble point 
characterized by d-isotopy. We remark that additional free- 
dom in defining Tr can be gained by exploiting the the ambi- 
guity of whether a bigon should be considered a loop or not, 
as discussed in L26tl . In particular, this allows one to make the 
diagonal elements of Tr equal. 

This construction shows how an extended Hubbard model 
with an additional ring exchange term (or the equivalent 
Quantum Dimer Model) can be tuned to the d-isotopy state(s). 
As discussed earlier, they satisfy two of the three conditions 
which define a class of stable, gapped topological phases 
which are centered about the special values d = 2 cos(7r/fc + 
2). The next step is to understand how perturbations can push 
the system (by implementing the JW projectors) into these 
phases. Our simplest candidate for a "universal quantum com- 
puter" is associated with d = (1 + -\/5) /2. 
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